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In non-central heavy ion collisions, identical two particle 
Hanbury-Brown/Twiss (HBT) correlations C(K, q) depend 
on the azimuthal direction of the pair momentum K. We 
investigate the consequences for a harmonic analysis of the 
corresponding HBT radius parameters ii^^-. Our discussion 
includes both, a model-independent analysis of these param- 
Leters in the Gaussian approximation, and the study of a class 
hydrodynamical models which mimic essential geometrical 
^j^d dynamical properties of peripheral heavy ion collisions. 
A|lsn we discuss the additional geometrical and dynamical in- 
f(| rmation contained in the harmonic coefficients of R^^. The 
fading contribution of their first and second harmonics are 
^^bu,nd to satisfy simple constraints. This allows for a min- 
Q^al, azimuthally sensitive parametrization of all first and 
^^cond harmonic coefficients in terms of only two additional 
fit parameters. We determine to what extent these param- 
can be extracted from experimental data despite finite 
Silultiplicity fluctuations and the resulting uncertainty in the 
^*«©construction of the reaction plane. 
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^ . I. INTRODUCTION 

' , The goal of the current and future experimental heavy 
Ciin programs at CERN and BNL is to test the equih- 
^ation properties of hadronic matter at energy densities 
ir^here quarks and gluons are the relevant physical de- 
J^ees of freedom. Anisotropic (directed) transverse flow 
kk an important observable for this program, since both 
<n]vdrodvnamic and thermodynamic behaviour is based 
CSi equilibration processes between local degrees of free- 
dom. Hydrodynamical flow effects result from pressure 
gradients which due to compression of the hadronic mat- 
ter build up during the collision process. Their strength 
depends on the equation of state of the hot matter and 
provides insight into the collision dynamics. Moreover, 
hadronic observables mainly test the final collision stage 
and a back extrapolation is needed to extract from them 
information about the hot and dense earlier stages. For 
this to work, collective and random ('thermal') motion 
in the collision region have to be distinguished properly. 
Hence, concepts of different types of collective flow play a 
central role in understanding the dynamics of heavy ion 
collisions. 

Directed anisotropic flow was observed in both 
AGS 1^,^ and SPS experiments, as well as at lower 
BEVALAC/SIS energies §. Directivity |,|, 2- and 3- 



dimensional sphericity , or the so-called deformation 
parameter Rp ^ are typical variables used in its charac- 
terization. With minor differences, all of them are sen- 
sitive to azimuthal anisotropics in the triple-differential 
particle distributions. The most complete experimentally 
feasible parametrization is obtained in a Fourier expan- 
sion in the azimuthal angle for different values of rapidity 
and transverse momentum, P,p],PpO) 
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Here, the azimuthal angles ip}^ allow for the determina- 
tion of the reaction plane and the harmonic coefficients 
Vn characterize the size of the total vector sum of trans- 
verse momenta {n — 1), the approximate elliptic shape 
of the azimuthal distribution {n = 2) and higher order 
triangle-type (n = 3), rectangle-type (n = 4) etc. defor- 
mations. 

In Eq. (1.1), we have expressed the one-particle dis- 



tribution in terms of the emission function S{x,p) which 
specifies the collision region at freeze-out. S{x,p) is a 
Wigner distribution and denotes the phase space prob- 
ability that a particle of four momentum p is emitted 
from a space-time point x in the collision region. Fea- 
tures of collective dynamics as e.g. directed flow are 
encoded in the source function S{x,p) as x-p position- 
momentum correlations. Observables extracted from the 



one-particle distributions (Id) are not sensitive to the 
space-time characterisitics (and a fortiori to x-p correla- 
tions) of the source. The question arises to what extent 
observables which are sensitive to x-p correlations can 
support and refine the picture obtained via the analysis of 
(pTl). This motivates an azimuthally sensitive Hanbury- 
Brown Twiss analysis of 2-particle correlation functions, 
which is the main focus of the present work. 

Identical two-particle correlations C(K,q), here writ- 
ten in terms of the average K = \{pi + P2) and rela- 
tive q = Pi — P2 pair momentum, are the only known 
observables sensitive to space-time characteristics of the 
source. Their space-time interpretation is based on the 
result |ll|-|l|l 



C(K,q) = l + |(e*i-("-''^*')|. 



/ d^xS{x, K) ' 
where we have used the on-shell condition q ■ K 



(1.2) 



1 



to substitute the temporal component in the 4- 
dimensional Fourier transform, (3 = 'K./Kq. According 
to (1.2), determining if-dependent geometrical and dy- 
namical source infomation reduces to a Fourier inver- 
sion problem (which due to the on-shell constraint how- 
ever does not have a unique solution). In the standard 
analysis of ( |1.2| ), one assumes an azimuthally symmetric 
collision and characterizes C(K, q) with four Gaussian 
Hanbury-Brown/Twiss HBT radius parameters which 
depend on |Kx| and the longitudinal pair rapidity Y 
only p^ . In contrast, the anisotropic case requires six 
HBT radius parameters Rij which in addition depend on 
the azimuthal angle $ of |jl^. Previous discussions 
of the azimuthal dependence of C(K,q) were based on 
event generator studies |17 1^ for finite impact parame- 



ter collisions, or exploited the Lorentz invariance of the 
correlator to derive azimuthally dependent HBT radius 
parameters pj|] . The present work is complementary to 
these and starts in analogy to ( |l.]| ) from expanding the 
angular dependence of C(K, q) in a harmonic scries, 

C(K,q) = l + exp {Y,Rlm}j , 

oo 

+ 2^i?^.„2(if^,r) cosn$ 

n=l 

oo 

+ 2Y,R!jJ{K^,Y)iimn<i>. (1.3) 

ri=l 

Here, the components i,j are given in the "out-side- 
long" os/-system where the relative pair momentum has 
a transverse out component qo parallel to the pair mo- 
mentum Ki^, a longitudinal long component qi along the 
beam and a remaining side component. To discuss the 
azimuthal ^-dependence of the two-particle correlator 
( |l.3| ), we frequently use an impact-parameter fixed sys- 
tem. In this system, the direction of the impact parame- 
ter b specifies the x-axis, the z-axis is along the beam and 
the 2/-axis is perpendicular to the reaction plane spanned 
by X and z. Accordingly, the total angular momentum L 
of the system, with 

= eijk{{xj Pk)) = eijk J J d*xxj pk S{x,p) , 

(1.4) 

points along the y direction. 

One central theme of the following is to determine 
those harmonic coefficients Rij^n^: Rtj.n^, whose contri- 
butions are not negligible. We shall find that there are 
very few independent ones. This makes a comparisons 
with experimental data feasible. Also, we aim at un- 
derstanding which geometrical and dynamical informa- 
tion about the particle emitting source is contained in 



the harmonic coefficients. In section |l[ we attack both 
problems by deriving model-independent expressions for 
the HBT radius parameters. These allow for the calcula- 
tion of HBT radii from ^-dependent space-time variances 
{XfiXi) of arbitrary model emission functions S{x,K). 
Investigating the ^-dependence of (Xf^x^) leads then to 
relations between the harmonic coefficients in (1.3). In 
the more detailed model-independent discussion in sec- 
tion III and the subsequent model study in section IV, 



we restrict our investigation to mid-rapidity and to sym- 
metric collision systems. Due to the reflection symmetry 
with respect to the y-z-plane, all odd harmonic coeffi- 
cients vanish in this case, and this considerably simpli- 
fies the discussion. In section 0, we extend this analysis 
to the fragmentation regions. Again, we derive simple 
relations between the non-vanishing first harmonic coef- 
ficients, and we illustrate our findings quantitatively in a 
subsequent model study. The discussion in sections 
implicitly assumes that the orientation ipn of the reac- 
tion plane is known. It hence neglects finite multiplicity 
fluctuations which introduce in practice a signiflcant un- 
certainty in determining ipn. In section VI we investigate 
to what extent information about the anisotropy of the 
correlator C(K,q) can be obtained despite these statis- 
tical constraints. The main results are then summarized 
in the Conclusion. 



II. AZIMUTHAL DEPENDENCE OF CARTESIAN 
HBT RADIUS PARAMETERS 

Here, we derive model-independ ent expressions for the 
Cartesian HBT radius parameters (|l.3| ) in terms of space- 
time variances [|l^,|l^ of the emission function S{x,K). 
These radius parameters depend in general on the az- 
imuthal orientation of the pair momentum K which we 
define with respect to the direction of the impact param- 
eter b, 



$ = Z(K_L,fo). 



(2.1) 



They can be calculated as second derivatives of the cor- 
relator C(K, q) with respect to the relative momentum 
components i,j = o,s,l in the os^-system. In what fol- 
lows, we consider the emission function S{x, K) to be 
given in the impact parameter fixed coordinate system. 
Then, to express the HBT radii in terms of space-time 
variances, one has to rotate the coordinate system by the 
angle <&, 



RUK) 



(3±\ fx cos ^ + y sin $ 

,X'$i= -isin<I> + ycos* 

a2c(q, K) 



dqt dqj 



= ([(P$.i), - (I?$/3),t][(I?$5)j - iV^I3)jt\) . 



(2.2a) 



(2.2b) 
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Here, = — (x^), and all coordinates cc, y and z 
are given in the impact-parameter fixed system. The 
space-time variances specify the curvature of the cor- 
relator at q = and coincide with the experimentally 
determined half widths of C(K,q) for Gaussian shapes 
only j2^,^l[]. Deviations from a Gaussian can be charac- 
terized by more refined methods |2l| . The present inves- 
tigation is restricted to correlators of sufficiently Gaus- 
sian shape and makes no attempt to quantify (possibly <&- 
dependent) non-Gaussian deviations. In the analysis of 
azimuthally symmetric HBT correlation radii, this Gaus- 
sian approximation has led to a qualitative and quantita- 
tive understanding of the Xj_-dependence of correlation 
functions jl5| . This motivates us to adopt the same start- 
ing point for an azimuthally sensitive analysis. The six 
^-dependent HBT radius parameters ( ^.2b| ) read 



Rl{K^, $, Y) ^ (i^) sin^ $ + {f) cos^ $ - {xy) sin 2$ , 
Rl{K^,<i>,Y) = (i2)cos2$ + (y2)sin2$ + /32 (P) 

-2/3±{ix) cos$ - 2P±{iy) sin(f> -f- (xy) sin 2$ , 
Rl,{K^,<i>,Y) = (iy) cos2$-f isin2$((y2) _ (^2^) 

+(3±_ (ix) sin <f> - (iy) cos $ , 
Rf{K^,^,Y) = {{S-f3ii)^), 

Rli{K^,'l>,Y) = ((z-/3if)(icos$-Hysin$-/3_Li)) , 
i?2^(ii'i,$,r) = ((z-/3;f)(ycos$-isin$)). (2.3) 

These equations separate the explicit ^-dependence of 
the HBT-radii (which is a consequence of the changing di- 
rection of the pair momentum K with respect to the reac- 
tion plane) from the implicit ^-dependence of the spatio- 
temporal widths {Xfj^Xy) (which reflects a ^-dependent 
change of the shape of the effective emission region). 
In general, both implicit and explicit ^-dependence will 
show up in the harmonic coefficients 



^J^l^f-.Rl cos(m<I))d<I>, 

s 2 
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^ Rl sin(m$) d$ , 



(2.4a) 
(2.4b) 



which determine the complete Gaussian parametrization 
(0). 

The following disc ussion of the ^-dependent HBT ra- 
dius parameters (2.3) is focussed mainly on the transverse 



parameters R^^, R^, R^^. Their harmonic coefficients 
depend on the ^-dependence of the transverse spatial 
widths 



global geometry of the collision region. For non-trivial 
x-^ correlations, this simple relation between the emis- 
sion region and the global geometry breaks down. In 
sections [II and 0, we find observable combinations of 
first and second harmonic coefficients which are sensitive 
to this difference. More generally, we classify the possible 



^-dependences of (2.5) and discuss their implications for 
the harmonic analysis of HBT radius parameters. 




FIG. 1. Schematic picture of a transverse cut through the 
collision region. Shaded is the effective emission region con- 
tributing to the correlator for pair momentum K. Two differ- 
ent scenarios are shown: in (a) the anisotropy is determined 
by the global geometry of the collision, in (b) it is significantly 
influenced by the collective dynamics of the source. Both sce- 
narios can be distinguished experimentally, see the text. 

We conclude this section by shortly commenting on 
the ^-dependence of Rf, Rli and i?^^. The longitudi- 
nal radius parameter Rf shows no explicit ^-dependence 
and coincides formally with the expression for the az- 
imuthally symmetric case [Q. The radii Rfi and R'^i 
contain explicit ^-dependent terms proportional to (tx) 
or (ty). These characterize asymmetries of the particle 
emission probability around the point of highest emissiv- 
ity and they vanish for models with Gaussian emission 
functions S{x, k). In the light of this, we consider in what 
follows all harmonic cocfffcicnts to > 1 of the HBT radii 



H ^ ^^oii ^si to be negligible. 



(2.6a) 
(2.6b) 
(2.6c) 



This is a reasonable but model-dependent assumption. 
For all models studied below, we have checked ( |2.6| ) nu- 
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merically, see section IVB 



(XiXj 



(2.5) 



i, j being components in the transverse plane. Any <&- 
dependence of is a consequence of non-trivial x-^ 
correlations and a fortiori of position-momentum corre- 
lations in the source. To illustrate this point, we have 
sketched in Fig. |l] two simplified scenarios. If there are 
no x-^-correlations, then the transverse shape of the ef- 
fective emission region is ^-independent and reflects the 



III. AZIMUTHAL ANALYSIS OF THE 
MID-RAPIDITY REGION 

In this section, we discuss the azimuthal analysis of 1- 
and 2-particle spectra for symmetric collision systems like 
Pb-Pb or Au-Au at mid-rapidity. The important simplifl- 
cation at mid-rapidity is that all observables are invariant 
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under $ — > $ + 7r, a 180°- rotation in the transverse plane. 
All odd harmonic coefhcients vanish. This is different for 
the fragmentation region or for non-symmetric collision 
systems where the only remaining symmetry is that with 
respect to the reaction plane. The arising complications 
are discussed in section 0. Here, we first discuss a sce- 
nario, for which the space-time variances (|]^) do not de- 
pend on the azimuthal direction of K . Then, we turn 
to the discussion of an arbitrary ^-dependence of T^. 



A. The elliptic approximation 



We start by considering an elliptic approximation of 
the transverse spatial widths . This toy example will 
be useful in the sequel for comparisons with the general 
case. It provides a simple picture for the consequences 
of a purely geometrical scenario, and it can account for 
some of the main features of the harmonic coefficients, 
calculated in the model study in section |rv| . 

In the elliptic approximation, we characterize the ten- 
sor T"*- by its principal axes (eigenvectors) 51, 52, and the 
corresponding eigenvalues gi and 52- If the orientation 
of the reaction plane and gi differs by an angle then 
takes a simple form in the impact-parameter fixed 



system. 



1 -|- ttT cos 2iy9 —aTsmlLp 
—ut sin 2lp 1 — aT cos 2ip 
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(3.1a) 
(3.1b) 



This is somewhat analogous to the expression of the 
transverse sphericity tensor in It provides a con- 

venient parametrization of the three independent expec- 
tation values (xiXj) in terms of the average size g of the 
homogeneity region, its transverse spatial anisotropy ax 
and the orientation tp of its principal axis gi with respect 
to the reaction plane. In case of an azimuthally symmet- 
ric collision with vanishing impact parameter, tp — 0, the 
cross term vanishes, {xgio) — 0, and the spatial asym- 
metry ax in ( 3.1a ), ( [3.1b| ) describes the difference be- 
tween the spatial widths in the out- and sirfe-directions. 
In terms of these parameters, the 'transverse' HBT radii 
Ro, Rs and Ros read 



(3.2a) 



i?2(A'i, $, y) = 5(1 - aT cos2((^ -I- $)) , 
i?2(X^,$,y) =g(l-HaTC0s2((^ + $)) 

+Pl{P)+Co, (3.2b) 

i?2,(X±,$,y) = -gaTsin2((^ + $) + Cos, (3.2c) 

Co = -2/3±(t5) cos $ - 2f-i^{ty) sin$ , (3. 2d) 
Cos = +/3± (fe) sin $ - /3i {iy) i 



> cos $ . 



(3.2e) 



The two correction terms Co, Cos contain widths {ixi) 
which are linear in t and measure asymmetries of the 
source around the point of highest emissivity. At mid- 
rapidity, all observables are invariant under $ $ -|- tt. 



and {xi) {(^ + tt) = — (if)($). The correction terms hence 
do not contribute to the first harmonic coefficients 
but to the second ones. Moreover, they vanish for emis- 
sion functions S{x, K) which are Gaussian in the spatial 
components. Their contribution to i?^, R^^ is neglected 
in the remainder of this section, and the validity of this 
approximation is checked in the numerical model study 



of section IV 



The main assumption of the elliptic approximation 



which we avoid in section III B , is that the parameters 



g and ax do not depend on the azimuthal orientation of 
K^. The homogeneity region is determined entirely by 
the geometry, cf. Fig. |lfL. The eigenvectors gi, g2 lie par- 
allel and orthogonal to the reaction plane and the angle 
if should vanish; it hence accounts only for the statis- 
tical uncertainty in determining the reaction plane, see 
Comp arin g the expressions ( |3.2| ) with the 



section 



VI. 



Fourier expansion (O) of the HBT radius parameters, 
we find for — Q 



Ro/^9 + mn 



R. 



0,2 



- -Rs,2 - 



Rs,0^ 
2 



■os,2 



-,axg- 



(3.3a) 
(3.3b) 



The other second and higher order Fourier coefficients 
vanish. According to (3.3), the ansatz ( |l.3| ) for the cor- 
relator contains redundant information. In the present 
case, neither the shape of nor the size of the homo- 
geneity region depend on $ while the anisotropy terms 
of the HBT radius parameters are sufficiently many to 
encode for such a non-trivial ^-dependence. 



B. The general case 

To discuss an arbitrary ^-dependence of the transverse 
spatial widths, we start from a complete Fourier expan- 
sion of Tjj^ = {xiXj), respectively of its linear combina- 
tions 

A = i£!>±il!l = J2n=0 i^n COS 71$ + siu n<i>) , 
B = ii!>^ = J2^^^ (S„ cos 77$ + sin 77$) , 



C = (xy) = ^ (C'j cos 77$ + Cn sin 77$) 



(3.4) 



ri=0 



Here, the coefficients A„, A'^, etc. are functions of the 
longitudinal pair rapidity Y and the transverse pair mo- 
mentum |Kx I , only. The symmetry of the collision region 
with respect to the reaction plane implies that the terms 
(i^), (y^) are invariant under $ — $, while the term 



{xy) changes its sign. Calculating the coefficients in (3.4) 
via Fourier transform, it is easy to check that the primed 
ones vanish. 



A' — B' — C' — . 



(3.5) 
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The general ^-dependent HBT radius parameters (2.3) 
read 



-Roc 



^-Bcos2$-Csin2$, 

A + Bcos2^ + C sin 2$ + Pi (P) , 

-Bsin2$ + Ccos2$, (3.6) 



where we have dropped the small correction terms Co, 
Cos- From these HBT radius parameters (3.6) one can 
calculate the harmonic coefficients (2.4). Especially, we 
obtain the zeroth harmonics 



iB2 + + I3l {P 



(3.7) 



and the second harmonic coefficients 



DC ^ 

Rs.2 - 



R 



OS, 2 



^Bq + jBi + jd , 



(3.8) 



For the case of ^-independent spatial widths T-^, only 
the terms with index survive on the right hand side, 
and ( |3.8D coincides with (3.3b) for Aq = g and Bq = 
OtQ- If deviations from ( 3.3b| ) are observed, this is an 
unambiguous sign for source gradients leading to a <i>- 
dependence of T^. From the absence of such deviations, 
however, one cannot conclude that there are no source 
gradients. In the following model study we shall find 
that even in the presence of sizeable source gradients, 
s uch d eviations can be small. Then, the use of equations 
( |3.3b ) resides in reducing the number of fit parameters 



in an azimuthal HBT-analysis, see section VI 



Let us finally anticipate that for the models studied 
below, the fourth harmonic coefficients A, B and C are 



ne gligib le. It follows from (3^) that then the deviations 
of ( 3.3b| ) are essentially determined by the term 1^2 only. 



I.e., 



r: 



■0,2 



2R: 



■OS. 2 



R 



S.2 



(3.9) 



IV. A MODEL CALCULATION 

We introduce now a simple hydrodynamical model for 
the emission function of a heavy ion collision which in- 
cludes anisotropy effects. For this model, we calculate 
both the one- and two-particle spectra, illustrating the 
main points of the above model-independent discussion. 

In the central rapidity region of a peripheral collisions, 
the initial distribution of the highly excited nuclear mat- 
ter is given by the intersection of the nuclear spheres. 



The largest pressure gradient developping from such ini- 
tial conditions is expected to be aligned with the impact 
parameter 5, cf. 18|. To mimic this scenario, we con- 
sider the class of model emission functions 



S{x, K) — Tom± cosh(?7 — y) exp ^ 



H{x), 
(4.1) 



whose azimuthally symmetric versions have been dis- 
cussed extensively in the literature. For a review, cf. 
These models assume the emission of particles from a 
thermalized system with collective four- velocity m^, con- 
fined in a space-time volume determined by H{x). The 
factor P ■ n(x) — TQm± cosh(77 — y) specifies a simple hy- 
perbolic freeze-out hypersurface. We introduce an az- 
imuthal anisotropy in the source both via an elliptic 
shape of the geometrical emission region. 



H{x) = exp 



{T-Tpf _ {1] - r/o) x'^ y"^ 
2(Ar)2 2(Ar;)2 2pl 2pl 



(4.2) 



and via an azimuthally asymmetric flow pattern u^{x) 
which is properly normalized, u^u^ — 1, 



Uf^{x) = {ui cosh r],Ux,Uy,ui sinh 77) 



— JL 



±, ui = Jl + ui + uf^ 



(4.3) 



In the longitudinal direction, we choose a boost invari- 
ant flow pattern satisfying Bjorken scaling |22|, i.e., the 
main energy flow is along the beam axis. Instead of the 
variables px, Py, Xx-, Ay, we use in what follows the trans- 
verse size R and its spatial anisotropy e^, as well as the 
transverse flow strength 77/ and the corresponding flow 
anisotropy e / , 



= RVT^s 



Py = R\/l + Es 



(4.4a) 
(4.4b) 



We have chosen the principal axes of the transverse space- 
time distribution aligned with those of the azimuthal mo- 
mentum distribution, since the dynamical evolution of 
the collision region cannot break the reflection symmetry 
of the system with respect to the reaction plane. The spa- 
tial anisotropy takes values in the range — 1 < es < 1, 
i.e., for e > 0, the collision region is longer in the direc- 
tion perpendicular to b, px < Py For the flow anisotropy 
£/, we allow for —1 < e/ < 1; the major flow component 
lies in the reaction plane if £/ is positive. All numerical 
calculations are done for the input parameters T = 150 
MeV, m = = 139 MeV, tq = 5 fm/c. At = 1 fm/c, 
A77 = 1.22 and R — 5 fm. We study the dependence of 
the 1- and 2-particle spectra on the size of the transverse 
flow, and the spatial and dynamical e/ anisotropies. 
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A. Harmonic analysis of azimuthal particle 
distributions 

The harmonic coefficients w„ of the triple-differential 



one-particle spectrum (1.1) are given in terms of the 
Fourier transforms 



dN 



cos{nq 
sm{n(j. 



r2TT 
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Vn cos(nV'^-') , 



(4.5) 



(4.6) 



The zeroth harmonic vq denotes the azimuthally aver- 
aged double differential particle distribution, all other 
anisotropy parameters Vn characterize azimuthal asym- 
metries in the momentum distribution. 

In the absence of transverse flow, rjf = 0, the model 
(4.1) does not contain source gradients in the transverse 
plane. Irrespective of whether the transverse geometry 
of the source is radially symmetric (e^ = 0) or not, par- 
ticle emission is isotropic in the transverse plane. All 
harmonic coefficients w„, n > 1, vanish. 

In the presence of transverse flow , the reflection sym- 
metry of the emission function (4.1) in the direction of 
the impact parameter implies a (f) (f) + n symme- 
try of the particle spectra. All odd harmonic coeffi- 
cients Vn vanish. The lowest non-vanishing anisotropy 
parameter is V2. This parameter is positive if the spec- 
trum N / (pt dpt dy d(f)) shows a maximum in the reac- 
tion plane, it is negative for the opposite case. Higher 
fourth (6th, etc.) harmonic coefficients are found to be 
much smaller. They are not discussed further. 

The transverse pair momentum dependence of the 
corresponding normalized quantity v^jv^ is depicted in 
Fig. ^for different physical scenarios. Irrespective of the 
model parameters f//, and e/, the coefficient V2/vq van- 
ishes at K±^ = and is growing monotonously with the 
transverse pair momentum. This is a direct consequence 
of the Lorentz-invariant Boltzmann term in ( [l.lj ) which 
encodes the assumption of local thermal equilibrium. 



K^u^ = ^ cosh(7, -Y)ui- - ^ 



(4.7) 

The terms which couple the flow components Ux, Uy lin- 
ear to the transverse momentum components are the only 
ones in the model emission function (4.1) which can in- 
troduce a ^-dependence. In the limit K± 0, these 
^f-dependent terms vanish, the emission probabilities in 
different azimuthal directions become equal, and 



lim Vn{K^) = 



for all n > 1. 



(4.8) 



More explicit expressions for the parameters Vn of this 
model can be obtained in a saddle-point approximation, 



details of which are presented in Appendix In this 
approximation, we flnd that V2 cx for small values of 
K±. For small transverse flow, the leading dependence 
on the anisotropy parameters and e/ is given by 



2(e/-e.) 



(4.9) 
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FIG. 2. The normalized second harmonic coefficients V2/V0 
of the single particle spectrum for the model ( [l.l[ ) . The diflter- 
ent plots show results for scenarios with diflterent transverse 
flow strengths rjf and geometrical anisotropies eg. The lines 
denote different transverse flow anisotropies: e/ — 0.8 (solid 
line), e/ = 0.4 (dashed line), e/ = (dotted line), e/ = —0.4 
(dash-dotted line) and e/ — —0.8 (thin solid line). Positive 
values of e/ correspond to the major flow direction lying in 
the reaction plane. 



Equation (4.9) describes correctly the main features of 
the numerical results in Fig. ||. The coefficient V2{K±) 
vanishes for = eg, and its sign coincides with that of 
tf — €s. In the case of an azimuthally symmetric particle 
emission region, — Q, this behaviour is entirely due to 
the transverse flow anisotropy ef. For positive values of 
the major transverse flow component lies in the reac- 
tion plane and V2 is positive. Negative values of mimic 
a squeeze-out scenario where the major flow component 
lies orthogonal to the reaction plane. This explains the 
e/-dependence oiv2, shown in Fig. ^. On the other hand, 
if we choose for a given flow pattern {rjf, ef flxed) more 
and more elongated transverse geometries, then V2 de- 
creases. The reason is that increasing ts results in more 
emission points with a large Uy and small Ux flow com- 
ponent. It thus mimics a larger squeeze-out component 
of the transverse flow orthogonal to the reaction plane. 
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In the simple model discussed here, increasing eg and 
decreasing e/ hence affects the azimuthal particle dis- 
tributions similarly. Spatial and dynamical information 
cannot be disentangled completely on the basis of single- 
particle spectra. For azimuthally symmetric scenarios, 
this is well-known [E3|. 



harmonics. Higher harmonics are found to be very small. 
Rf shows no ^-dependence, the only non-vanishing sec- 
ond harmonics are R^i ^ and i?^; ^ . From Fig. ^ we con- 
clude that these higher harmonics are negligible. This 
illustrates the arguments leading to (2.6). It suggests 
to restrict an azimuthally sensitive HBT analysis to the 
'transverse' HBT radius parameters i?^, i?^ and i?^^. 



B. Harmonic analysis of HBT radius parameters 



Here, we check first that the ^-dependence of the radii 
i?^; , and i?^; is negligible. At mid-rapidity, Y — (ji — 



0, the emission function (4.1) is symmetric with respect 
to z ^ 
vanish, 



-z, and the HBT radius parameters i?^; and R\ 



Rli,Rl^{{S-Pii)), 

{{i - Pii))\Y=o = foraUKj.. 



(4.10) 



rij = 0.3, 8, = 0,63=0.1 
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efficients of the non-transverse HBT radius parameters, R, 
(solid line), Roifi'^ (dashed line), 7?^; 2^ (dotted fine) and 
Rli 2^ (dash-dotted line) for forward rapidity Y = 1. The sec- 
ond harmonic coefficients R'01^2^ and ii^; 2^ &re proportional 
to {x) and (y), which are generically small correction terms, 
see the discussion following ( 4.1C ) . 
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FIG. 4. Zeroth and second harmonic coefficients for the 
transverse HBT radius parameters of the model ( [l.l| ) with- 
out (77/ = 0) and with (77/ = 0.3) transverse ffow, and for 
different sizes of the spatial anisotropy e^. Thin and thick 
solid lines denote the zeroth coefficients i?^ q s-nd R'ifi re- 
spectively. The plot shows all nonvanishing second harmonic 
coefficients, Ro.2'^ (dash-dotted line), R's^2^ (dashed li ne) a nd 
^03,2^ (dotted line). For rjf = 0, 7?^, 



Rls,2 , see (3.3b) 



To study their K^-dependence, we hence choose the 
forward rapidity Y = \. Numerical results are pre- 
sented in Fig. ^ The zeroth harmonic coefficients show 
the i^_L-dependence expected from azimuthally symmet- 
ric model studies Especially, the longitudinal ra- 
dius parameter R^^^ has a steep ifj^-slope which refiects 
the strong longitudinal expansion of the source. Also, 
the out-longitudinal cross term i?^; q shows the typical 
iC_L-dependence known from studies of azimuthally sym- 
metric models. It takes significant non-zero values and 
vanishes at K±^ = where R^^ = i?^;. The parameter 
i?^; Q vanishes for all Ki_ . For an azimuthally symmet- 
ric emission region, this would be a consequence of the 
qs reflection symmetry of the correlator. Contri- 
butions breaking this symmetry introduce automatically 
a ^-dependence and hence do not show up in the zeroth 



For the transverse HBT radii of the model (4.1), sug- 
gestive analytical expressions can be obtained in the ap- 
proximation w 1 -fju^ + 5"^- Deferring all technical 
details to Appendix we merely state that in this ap- 
proximation, the (f>-dependence of T-*- is lost. The zeroth 
and second harmonic coefficients of all transverse HBT 
radii can be written in terms of an average size g and 

Eq. 



a spatial anisotropy ar, introduced in section |n 



gar 



^R2 (,^+,^I^^2(i_^2)^ 



l 2^ ?7/ H 



4) 4 
—"If 



(4.11a) 
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i?2 1 



9 = 



1 



f, mj_(l-eae/) ^2 



T2 ?7/ 



(4.11b) 



In the hmit of vanishing transverse flow, t^/ ^ 0, these 
expressions become exact. The value of g which deter- 
mines the zeroth harmonic of is just given by the 
transverse geometrical radius R^. The anisotropy gar 
which specifies the second harmonic coefficients, is pro- 
portional to e^. This is the case shown in Fig. ^. For 
77/ — 0, the model emission function has no intrinsic 
position-momentum correlations, and the different har- 
monic coefficients satisfy the relations (3.3). The com- 



ponents R'g2 ^''^d 2 coincide and they differ from 
i?Q 2^ by an overall sign only. 



In the presence of realistic transverse flow strengths 
the approximation ui \ -\- -|- looses its va- 
lidity. Numerical calculations are needed to make precise 



quantitative statements, but the expressions (4.11) still 
describe essential qualitative features. In the limiting 
case of an azimuthally symmetric collision region with 
a finite transverse flow 77/, the anisotropy parameter 
vanishes and g goes to the well-known ||l^ lowest order 
expression for the side radius parameter, 



Rs'^-g = 



i?2 



(1 + ^,^) 



fore^^e/^O. (4.12) 



This describes the leading mj^-dependence of Rsfl" as a 
function of 77/: the slope of the side radius parameter is 
indicative of the transverse flow rjf. 

In the absence of dynamical anisotropies (e/ =0) the 



second harmonics R'i 



0,2 



R 



S.2 



and R 



OS. 2 given essen- 



tially by ^gar, are sensitive to the strength of the geo- 
metrical anisotropy es . For small values of 77^ , the leading 
Cg-dependence of ( 4.11a ) is linear and this is consistent 
with the scaling of the dashed and dash-dotted (e/ = 0) 
lines in Fig. ^. The if i^-slopes depicted in Fig. ^ are quali- 
tatively explained by ( 4.11a ). Also, we have investigated 
numerically the e/-dependence of the HBT radius pa- 
rameters for the model of section IV A. Here, we merely 
state that the main qualitative features of the numeri- 
cal results can be understood in terms of the analytical 
approximations ( 4.11a ), ( 4.11b ). 



It is a remarkable feature of the model (4.1) that the 



transverse spatial widths calculated in a saddle-point 
approximation do not show any ^-dependence. This can 
be traced back to the the ^-dependent terms in (4.7) be- 
ing linear in x and y, rather than e.g. quadratic. The 
slight difference R^ 2^ ^ R^s 2 found in the numerical 
calculation of Fig. |^ for r]f — 0.3 stems from the error 
made in the Gaussian approximation. To avoid drawing 
conclusions on the basis of a very model-dependent fea- 
ture like this almost complete cancelation of ^-dependent 
contributions for linear flow profiles, we have investigated 
the following two non-linear flow profiles as well, 



[square root profile], (4.13a) 



Ui ^ 77/-\/l ± £/ "^'j^l'^ , [quadratic profile]. (4.13b) 

Here the index i runs over x, y. Results for vanishing 
flow anisotropy e/ are shown in Fig. |^. In the limit 
— !■ 0, the dependence of on the azimuthal direc- 
tion of has to vanish and hence, all scenarios depicted 
in Figs. ^ and |^, confirm the purely geometrical relation 



O, R, 



?Q 2 — —Rs 2 — ~Ros,2 ■ For increasing val- 
ues of K±, we find deviations from this relation which 
are clearly more significant for the non-linear flow pro- 
files. These howe ver are very well accounted for by the 
modified relation ( 3^ ) amongst the second harmonic co- 



efficients, Rg 2 +2-Ros,2 = Rs,2 ■ The deviations from 
a purely geometrical scenario seen in Fig. ^ are hence not 
due to the correction terms Co, Cos, which are generically 
negligible, but to the leading ^-dependence of T-^ which 
enters (3.7) via the term A2. 




flow profile: square root 



10 - 



10- 




flow profile: quadratic 



10 



200 



400 



RI2 is still valid, see Eq. (3.9) 



600 800 
Kj^ (MeV) 

FIG. 5. Same as Fig. ^ but for a square root and a 
quadratic transverse flow profile, for which a Gaussian ap- 
proximation of the emission function is not ^-independent. 
In these cases, the difference R"^ 2^ ~ Ros,2^ is more signifi- 
cant, but i?^_2^ +2i?L,2^ 
and text below. 



For the sake of completeness, the zeroth harmonic co- 
efficients are presented in Fig. ^ as well. Details of their 
i^x-dependence can be understood by investigating the 
dependence of (x^) and (y^) on the pair momentum p^ . 

Let us sum up the results of our model study: All 
second harmonic coefficients of non-transverse HBT ra- 
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dius parameters are negligible. Amongst the transverse 
HBT radii, three second harmonic coefhcients are non- 
neghgible, namely R%^-i ^ -^^,2^; ^'^'^ ^os,-! ■ Their lead- 
ing contribution is determined by one parameter only, see 



( ^.3b ). The leading deviation from (3.3b) satisfies 



and is indicative of strong correlations in the source. 



V. AZIMUTHAL ANALYSIS OF THE 
FRAGMENTATION REGION 

For finite impact parameter, the collision system has 
a finite total angular momentum L, Li ~ Cijk{{xj pk)) . 
In heavy ion collisions, L and Ltot = \L\, is not di- 
rectly observable but has to be infered indirectly on the 
basis of impact parameter dependent measurable quan- 
tities like total particle multiplicities or transverse en- 
ergy. In general, angular momentum conservation is a 
constraint which can affect the shape of two-particle cor- 
relations p5|] . 

For the incoming nuclei, the angular momentum is en- 
tirely determined by Ltot = ((a^Pz )), but in the final state, 
the component ((z px )) in general carries part of Ltot ■ To 
illustrate the consequences, we consider a longitudinally 
expanding collision scenario, for which the space-time ra- 
pidity of the emission points ^7 = ^ log is linearly 
related to the momentum rapidity Y of the emitted par- 
ticles, rj = rji ■ Y. Using z = rsinhr/ = Tsinh(r/;y), 
one sees that a non- vanishing contribution {{zpx )) to the 
total angular momentum leads to a non-vanishing total 
vector sum of px in non-central rapidity bins, and hence 



to a non- vanishing Fourier coefhcient vi in (1.1). This ef- 
fect is more enhanced for larger rapidities Y, it vanishes 
at mid-rapidity Y = and it has opposite sign for back- 
ward rapidities. The size of the fraction {{zp^ )) of Ltot, 
depends on details of the dynamics. Note however, that 
((z Px)) 7^ does not imply necessarily the equilibration 
of microscopic degrees of freedom or source gradients. 

According to this simplified scenario, the 180°- 
symmetry of the transverse collision region is lost at non- 
central rapidities, and the odd harmonic coefRcients of 
the one- and two-particle spectra do not have to vanish 
any more. In this section, we investigate the implications 
for the first harmonics of the HBT radius parameters. 



A. Properties of first liarmonics 



To calculate odd harmonic coefficients of the trans- 



verse HBT radius parameters (2.3), we start again from 
the harmonic expansion (3.4) of T^. It is a remarkable 



property that in the radii Rjj ^ , Rfj ^ , even and odd 
harmonic coefficients An, A'n, etc. do not mix. For m 
even, only terms with n even appear, and for m odd, 
only terms with n odd. Especially, we find for the first 
harmonic coefficients 



^s,i^ = 5^1 ~ 3-^1 ~ iC"! - - , (5.1a) 

- i^i + \Bi + ICi + + IC3 , (5.1b) 

Ks/ = -3^1 - jCi + is, + IC3 , (5.1c) 

= K,,^ ^R^s/ = 0. (5.1d) 

It follows immediately that all first harmonics vanish in 
the absence of source gradients. 



r: 



= R 



0. 



[no source gradients]. (5.2) 



The physical reason is that for a source without gradients, 
the effective emission region has the same side and out 
extension irrespective of whether it is viewed under an 
angle <i> or an angle ^ + tt. 



For the case that all third harmonic coefficients in (5.1) 
are negligible, the three non- vanishing HBT radii in (5.1) 
satisfy 



pc 2 



2m 



(5.3) 



This equation is reminiscent of (3.9). There, however, 
the contribution Bq is typically an order of magnitude 
larger than the second harmonic contribution A2, and 
this allows for the further simplification (3.3b). Here, in 
contrast, all leading terms are first harmonics. The first 
harmonic coefficient of (i^) should be much larger than 
that of (y^) since asymmetries with respect to the beam 
axis will occur in the direction of the impact parameter 
only. In this limiting case which is relevant for the models 
studied below, we have Ai ~ Bi ^ Ci, and 



R' 



Rt 



Ri 



1 



(5.4) 



B. Model extensions for non-central rapidities 



The model emission function (4.1) investigated in sec- 
tion IV describes a collision with vanishing angular mo- 
mentum. Here, we introduce a simple extension which 
allows for finite angular momentum Ltot but coincides 



with the model (4.1) at mid-rapidity. Clearly, such ex- 
tensions are numerous. One can e.g. shift the transverse 
flow pattern in the x-direction, 



X + xy 

R 



(5.5) 



Since x multiplies the rapidity Y , the flow pattern shows 
a forward-backward anticorrelation in Y and will result in 
non- vanishing o dd c oefficients w„ and a finite total angu- 
lar momentum (|l.4| ). Alternatively, one can introduce a 
rapidity-dependent deformation of the transverse geom- 
etry, e.g. by i ntroducing a dependence of the Gaussian 
widths in (4.2) on the azimuthal direction tp, x = r cosip 



Px = (R + xycosif) Vl ~ 
Py = {R + xY cos (p) VI + 



(5.6) 
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For finite transverse flow, this again results in non- 
vanishing odd coefficients w„ and a finite Ltot ■ Both these 
modifications, ( ^.5[ ) and (5.6) reduce at mid-rapidity 
F = to the model studied in section [V. We have 



X = 2fm 



Y = 1 



X = 4 fm 



investigated the iir_L-dependence of their first and sec- 
ond harmonics and the (approximate) relations which 
they satisfy. The linear coupling of ( ^.5| ) on 'f'-dependcnt 
terms in (4.7) makes the x-dependence of the model (5.5) 
very weak. (The term K^u^^ does not introduce an 
additional position-momentum correlation, and the only 
X-dependence of the correlator stems from ui). Here, 



we present results for the model (5.6) only. To isolate 



the effect of the x-displacement, we have set the other 
anisotropy parameters to zero, eg = e/ = 0. All calcula- 
tions are done at forward rapidity, Y = 1. 

For vanishing transverse flow, rjf = 0, our model con- 
tains no source gradients. The zeroth harmonic coef- 
ficients show the expected behaviour: Rgfi is a K±- 
independent constant from which the out radius i?o,o 
differs by the factor P\{t^) only. Also, in the ab- 
sence of source gradients, all first harmonics vanish, 
and (5.2) holds. The second harmonics do not vanish: 
they are if_L-independent parameters determining the 
elliptic approximation of the transverse source geome- 
try. Also, they satisfy the geometrical relation ( ^.3b| ) 
as expected for sources without position-momentum cor- 
relations. Qualitatively, their behaviour is completely 
consistent with the case discussed for 77/ = in Fig. ^. 
Quantitatively, we find for the non- vanishing components 



R 



■y",2 



0.5 if X = 2 fm and \R. 



1.5 if X = 2. 

This is slightly larger than the values obtained for the 
case 77/ = 0.3 at K±^ = (see Fig. |^) and the minimal 
r7/-dependence can be traced back to the wj-dependent 
term in the Boltzmann factor (4.7). 

For finite transverse flow rjf — 0.3, numerical results 
are presented in Fig. ^. Due to source gradients, the az- 
imuthal eccentricity introduced via ( |5.6|) now shows up in 
the first harmonic coefficients. They vanish at vanishing 
transverse pair momentum. 



Jim R^^/{K^, 



y) = 0, (5.7) 



since the correlator is at K±^ = sensitive to the ge- 
ometry of the source only. More importantly, over the 
complete -fir_L-range, the first harmonics in the out, side 
and out-side directions satisfy up to a few percent the 
relation 3 : 1 : —1 as we have argued in deriving (5.4). 
The second harmonic coefficients satisf y for K±^ = the 
relation — 1 : 1 : 1 as suggested by ( 3.3b ). For finite 
K±, dynamically introduced deviations are found, but 



the modified relation (^, i?^ 2 + 2i?: 



scribes all results up to a few percent. 
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FIG. 6. Zeroth, first and second harmonic coefficients of 
the out (dash-dotted lines), side (dashed lines) an d ou t-side 
(dotted lines) HBT radius parameters for the model ( |4.l| ) with 
the modification (5.6) at forward rapidity Y = 1. To high 
accu racy, the first harmonic coefficients satisfy the relation 
(5.4) and the second harmonics the relation (BJJ). 



VI. FINITE EVENT STATISTICS AND FINITE 
MULTIPLICITY FLUCTUATIONS 

For the particle multiplicities obtained at the AGS and 
SPS, a determination of HBT radius parameters on an 
event-by-event basis is not possible. Finite event statis- 
tics limits the possibilities of a multidimensional HBT 
analysis for the now typical samples of the order of 10^ 
- 10^ events. To extract despite these statistical con- 
straints at least the major anisotropic HBT character- 
istics, it is clearly helpful to start from an azimuthal 
parametrization of C(K,q) in terms of a minimal set 
of fit parameters. Such a parametrization is discussed in 
section 



VIA 



In contrast to constraints from finite event statistics, 
the statistical uncertainties stemming from finite multi- 
plicity fluctuations cannot be overcome by investigating 
larger event samples. They constitute a fundamental lim- 
itation to any investigation of anisotropy measures, and 
they are particularly important in reconstructing the re- 



action plane. In section VI B we investigate to what ex- 



tent this affects the determination of anisotropy measures 
from two-particle correlators. 
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A. A minimal azimuthal parametrization 



B. Event samples with oriented reaction plane 



On the basis of our discussion in the sections || - 0, 
we propose as minimal parametrization of a two-particle 
correlator for peripheral collisions the Gaussian ansatz 

C^,(K,q) « 1 + A a^„(K,q) Ci(K,q,Vi?) 

xC2(K,q,Vfl) (6.1a) 
Csyni (K, q) = exp[-i?o^o^ ql - Rs,o'^ ql 

~Rl/qf ~2Roi/qoqi] (6.1b) 
Ci (K, q, tPr) = exp [ - ai (3 ql + q'i) cos($ - ?/'fl) 

+2aiqoqsSin{^-ijR)] (6.1c) 
C2(K,q,i/'i?,) = exp [~ a2{ql - g^)cos2($ - ipa) 

+2a2gogsSin2($- Vfl)] . (6.1d) 

The zeroth harmonic coefficients in this parametrization 
provide the complete Cartesian parametrization for the 
azimuthally symmetric case. Especially, the cross term 
Roi,o^ vanishes in the CMS at mid rapidity, and due to 
the symmetry arguments made in section the pa- 
rameter ai vanishes under these conditions as well. The 
parametrization for the first effective harmonic coefficient 
a I a nd th e seco nd one a2 is motivated by the relations 
(^) and ( |3.3b| ), i.e., it is based on setting 



a2{KuY) ^ Rl/ ^ -Rl/ ~ ~R 



os.l ' 
2 

OS. 2 • 



(6.2a) 
(6.2b) 



In the model studies in section |^ and 0, we have as- 
sumed that $ = corresponds to the pair momentum 
K lying in the reaction plane. Experimentally, this re- 
action plane is unknown a priori. The parameter ipR 
which determines its orientation has to be included in a 
comparison with experiment. A discussion of the statis- 
tical uncertainty in its determination and the implication 
for the extraction of the anisotropy parameters ai, a2 is 
given in section 



VIE 



In the model-independent analysis of the ^-dependent 
HBT radius parameters in sections 



III and 



IV, we have 



argued that corrections to these relations can be expected 
to be comparatively small on general grounds. Also, we 



have quantified deviations from (6.2) in the model studies 
of section IV and According to these studies, the most 
r eas onable non-minimal extension of the parametrization 
( |6t| ) is to implement the constraint i?^ 2^ +2i?; 



■OS, 2 



R1 2 instead of (3.3b), i.e., to replace (3.2b) by two pa- 
rameters for the second harmonics. 

Clearly, it would be preferable to fit to the most gen- 



eral parametrization (1.3), and to quantify corrections to 
these relations. However, as long as finite event statistics 
forces one to restrict the space of fit parameters as far 



as possible, the relations (6.2) provide in the light of our 
analysis the most reasonable set of constraints to adopt. 



The main problem in extracting the fit parameters ai, 
a2 from experimental data is, that due to finite event 
multiplicity, a multidimensional analysis of two-particle 
correlations is not possible on an event-by-event basis. 



Fits have to be done on event samples and the use of (6.1) 
presupposes that the different events are sampled with 
a fixed orientation of the reaction plane. The reaction 
plane can be reconstructed from an azimuthal analysis 
of the single particle distribution (IT), but its event-by- 
event determination is subject to significant statistical 
uncertainties. Since the parametrization (6T) depends 
on the angle these uncertainties affect the determi- 
nation of the HBT radius parameters as well. Here, we 
investigate quantitatively to what extent these uncertain- 
ties affect the determination of the anisotropy parameters 
ai, a2 from the experimental data. 

Our starting point is the assumption that the proba- 
bility distributions W{vi,tpii) of the first harmonic co- 
efficients of the one-particle spectrum around {vijiPr) is 
given by §§ 



1 



exp 



2viVi COSIpR 



2a2 



(6.3) 



We consider the case that the reaction plane is recon- 
structed from the first harmonic coefficients only, and we 
use the shorthand ipR = ipR^^- For ideas about improv- 
ing this reconstruction by taking higher order harmonics 
(wijV'i^'') i'^to account as well, we refer to ||]. Here and 
in what follows, we orient the most likely, 'true' direction 
along the x-axis, TpR = 0. The Gaussian distribution 
( |6^ ) can be used if event multiplicities N are sufficiently 
large to apply the central limit theorem. The variance 
then scales like , and the reaction plane is well defined 
for extremely large event multiplicities. 



lim W{vi,iPr) = S {iJjr) . 

N—*oo 



(6.4) 



For finite multiplicities N in the hundreds, however, the 
uncertainty in the eventwise determination of ipR cannot 
be neglected. The event sample with orien ted reaction 
plane should not be compared directly to (6T), but to 
an effective correlator which takes the probability distri- 
bution W of experimentally determined reaction plane 
orientations properly into account, 

C|«(K,q) = J v,dvid^PRW{vi,^PR)G^AK,q) . (6.5) 



It follows from the form of (6.3) that this effective cor- 
relator docs not depend on vi and cr separately, but is a 
function of ^ = vi/a only. The parameter ^ is a direct 
measure of the accuracy for the reaction plane orienta- 
tion [uHol . From the investigation of Voloshin and Zhang 
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(see Fig. 4 in |^), we conclude that a value of f « 2 cor- 
responds to an uncertainty of approximately 30°. This 
is the current experimental standard which will improve 
with the larger event multiplicities at RHIC and LHC. 

We investigate now in a simplified example to what 
extent finite anisotropies of the single events leave traces 
in an event sample constructed as outlined above. To 
this aim, we calculate for given anisotropies ai, a2 the 
effective correlator Cj*^ which is then fitted to the Gaus- 



sian ansatz ( |6.l| ). The anisotropy parameters (ai), (02) 
determined in this fit are then compared to the param- 
eters ai, a-2, of the single events one has started with. 
In Fig. the fitted averages (a^) are shown as function 
of the statistical uncertainty ^ in the reconstruction of 
the reaction plane from first harmonic coefficients. For a 
realistic uncertainty of 30° « 2), approximately 85 % 
of the 'true' first anisotropy parameters ol\ and approx- 
imately 55 % of the second anisotropy parameter ai is 
obtained. These ratios are independent of the size of a\ 
and 0L2- 

The correlator of an unoriented event sample is ob- 
tained from (3.1) by averaging C^^ over all orientations 



of the reaction plane. It is known that in fits to such 'un- 
oriented' correlators, HBT radius parameters receive ar- 
tificial contributions due to the averaging procedure . 
Since we average in C^J ^ over different reaction planes, 
such artificial contributions exist for C^*^ too. Especially, 
the zeroth harmonic coefficients of the HBT radius pa- 
rameters should show a ^-dependence. In the numerical 
analysis of the curvature of the correlator, we found this 
to be negligible. To understand the reason, we consider 
the case a2 = when the integral ( |6.5| ) can be evaluated 
analytically. The ?/'j?-dependent part reads 

Cf(K,q) = y" vxdvr dVflI^(«l,^i?JCi(K,q,7/^;j) 

idi^^^{~\{e^~e))h{Z). (6.6) 

The argument Z of the Bessel function depends on 
the relative pair momentum components go, Is and the 
anisotropy parameters a\ 

Z = {li)^ - 2ee {F, cos($) + Fs sin(a>)) 



-4 sin($) cos($) 



-2aiqoqs 



(6.7a) 
(6.7b) 



The limit ^ — > in (6.5) corresponds to an unoriented 
correlator. For this case, the argument Z is quadratic in 
the components of q and expanding Iq for small argu- 
ments, we find, (|A5|), 



d'lo(Z) 



dqi dqj 



= for ^ = 0. 



(6.8) 



q=0 



This illustrates for a simple example that the unweighted 
averaging over different event orientations discussed in 



|l9[ does affect the shape of the correlator in q, but not its 
curvature. It hence has to be determined by quantifying 
the deviations of C|*^ from a Gaussian shape [^. Here, 
we do not pursue this point further. 



~6 

E 

""5 




FIG. 7. The HBT anisotropy parameters (oi) as a func- 
tion of the parameter ^ which characterizes the event-by-event 
reconstruction uncertainty in the orientation of the reaction 
plan e. The parameters (qi), (0:2) are determined by fitting 
(6.1) to an event sample (6^) of correlators whose reaction 
planes are oriented alo ng the a;-axis according to the proba- 
bility distribution (6^). The value €, = 2 corresponds to a 
reconstruction uncertainty of approximately 30°. 



The calculation leading to Fig. ^ is based on simpli- 
fied assumptions. Especially, we have not considered 
eventwise fiuctuations in the size of the ai. Since both 
anisotropy fit parameter can take positive and negative 
values, such fiuctuations do not fake anisotropy signals. 
We hence conclude from the above that with a typical 
30° resolution of the reaction plane orientation, a de- 
termination of the anisotropy parameters ai, ai from 
experimental data is feasible. 

VII. CONCLUSION 

In the present work, we have studied the possibilities 
of an harmonic analysis of two-particle HBT radius pa- 
rameters. We have clarified to what extent the two fun- 
damental problems of such an analysis can be overcome: 

Firstly, the harmonic Fourier expansion of HBT ra- 
dius parameters introduces a plethora of additional fit 
parameters which make a comparison with experimen- 
tal data difficult. We have argued that under reasonable 
assumptions on the geometry and dynamics of the col- 
lision region, only very few of them are non-negligible. 
For the first harmonic coefficients, these are R'^i^ ^ , 
Ros,i^t and they can be described by one single parame- 
ter since their leading contributions scale like 3 : 1 : — 1 . 
Amongst the second harmonic coefficients, only R'^i^y 
Rg 2^, Ros -1 a-re non-negligible and their leading contri- 
butions scale like —1 : 1 : 1. Higher harmonic coeffi- 
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cients were found in all studies to be an order of mag- 
nitude smaller. Based on these observations, our main 



result is the parametrization (6.1) which describes the 
leading anisotropy of the two-particle correlator by two 
additional fit parameters only. If deviations from this 
parametrization turn out to be important, than this will 
provide an important constraint on further model stud- 
ies. A first non-minimal parametrization suggested by 
our studies quantifies deviations of the second harmonic 
coefficients from — 1 : 1 : 1 by one additional fit param- 
eter. 

The second fundamental problem in the determina- 
tion of anisotropy signals from C(K,q) is the statistical 
uncertainty in the eventwise reconstruction of the reac- 
tion plane. We have shown that for event samples with 
a typical 30° uncertainty in the eventwise orientation, 
the major part of the anisotropy measures ai, ai sur- 
vives. An analysis o f exp erimental data on the basis of 
the parametrization (6.1) seems hence feasible. 

We finally discuss which physical information can be 
extracted from the anisotropy parameters ai, oi2- From 
the discussion in sections |ll| - 0, we conclude that a 
non- vanishing first harmonic coefficient a\ automatically 
implies the existence of dynamical source gradients, see 
Eq. ( |5.2| ). In contrast, the leading contribution to a-i is 
determined by the geometry of the collision region and 
specifies essentially the elliptic shape of its transverse ex- 
tension. More explicitly, in the context of the hydrody- 



namical model of section [II, there is a tentative strategy 



to determine the geometrical and dynamical anisotropics 
Eg, e/ respectively. According to ( |4.9| ), experimental data 
on the harmonic coefhcient V2 restricts the allowed pa- 
rameter space in {es,ef) to a one-dimensional one. Ac- 



cording to (4.11a), the second harmonic coefhcients of the 



transverse HBT radii show a different dependence on 
and e f and can hence be used to constrain the remaining 
freedom. This illustrates that in non-central collisions, 
as in azimuthally symmetric ones, only a combination of 
1- and 2-particle spectra will allow to disentangle geo- 
metrical and dynamical information. 
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APPENDIX A: CALCULATION OF AZIMUTHAL 
PARTICLE DISTRIBUTIONS AND HBT RADIUS 
PARAMETERS 

In this appendix, we give details of how to calculate the 
harmonic coefficients of the single particle distributions 
( pTl] ) and the two-particle correlations (1.2) for the model 
in section fV, 



cosncj) / d'^x S{x, K) 



= Tom± J VriGniri) , 



(AI) 



Vt] ^ I T dr dr] coshirj - Y) e 2(Ar)2 2{A^y\ (^2) 



For n odd, all coefficients Vn vanish due to the 180°- 
symmetry in the transverse plane. In the approxima- 
tion Ul K 1 + ' -^"'^ 

quadratic in x and y, and 



2Uy, the Boltzmann term K^^u^ is 



Gniv) 



COS ncj) / dxdye 



T 



2P^ 



V 

2P! 



where we have used 

5x = 1/ f A + 4 1 . A 

\ Px 



2T2 



-I^{Z), (A3) 



-gy 



Kl 



= ^ cosh(77 - Y) , 

(K An 



2T2 



(A+f)(A+S) 



(A4) 



The modified Bessel function Iil in (A3) is obtained by 
doing the (/)-integral. To extract the leading dependence 
in Z, we expand for small arguments. 



Z COS 20 



E 

k=Q 



^ fc!r(f -f /c 



1 



(A5) 



For the coefficient n = 2, the leading Z-dependence in 
( A4) is linear. This implies (O). Also, for small K^^ 
A is approximately constant and Z oa . Since A de- 



pends on 77, the 77-integration in (AI) has to be done 
numerically. To obtain simple expressions for the main 
qualitative features, one may use 



(A6) 



This amounts to cosh(77 — y) « 1 and allows for the ap- 
proximate analytical calculation of HBT radius parame- 
ters. The latter are calculated via space-time variances 
which we express here in terms of the averages 
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{f{x,v))* = / dxdy f(x,y)e 



2pI 



2^ 



J Vrjjxfj.x^)* 



(A7) 



(A8) 



1„,2 



u?,, the 



Again, in the approximation ui w 1 , -^u.^ ^ ^^yi 
Boltzmann term K^u^_ is quadratic in x and y, which 



allows for an analytical calculation of the x- 
integration. We find {xy)^, = 0, and 



and y- 



R\l-es){l)^' 



l + ^7?2(l-e,)(l + e/) ' 
l + A7^){l + es){l-ef)- 



A- 



Py 



(A9) 



These averages are the building blocks for the space-time 
variance ( |A8|) which determine the HBT radius parame- 
ters. Since A depends on -q, the remaining r^-integration 
has to be done numerically. With the help of t he a pprox- 
imation (A£), however, the ry- integration in (A8) drops 
out and the analytical expressions ( 4.11 ) for g and gar 
can be obtained. Their validity is subject to the ap- 
proximations made above but they give a qualitatively 
correct, simple and intuitive description of the numerical 
results. Most remarkably, these approximate expressions 
do not depend on the azimuthal direction of K, cf. the 
discussion in section IV B. 
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